Experimental constraints on non-linearities induced by two-photon effects in elastic 

and inelastic Rosenbluth separations 
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The effects of two-photon exchange corrections, suggested to explain the difference between mea- 
surements of the proton elastic electromagnetic form factors using the polarization transfer and 
Rosenbluth techniques, have been studied in elastic and inelastic scattering data. Such corrections 
could introduce e-dependent non-linearities in inelastic Rosenbluth separations, where e is the vir- 
tual photon polarization parameter. It is concluded that such non-linear effects are consistent with 
zero for elastic, resonance, and deep-inelastic scattering for all and values measured. 

PACS numbers: 13.40.Gp, 25.30.Fj, 12.20.Fv 



I. INTRODUCTION 

Electron-proton (e-p) scattering has proven to be a 
powerful tool in the investigation of the structure of the 
nucleon. This interaction is typically described as the 
exchange between the electron and the proton of a single 
virtual photon. Because the electron is a point-like par- 
ticle, any structure observed in e-p scattering must be 
related to the target structure. Moreover, the relatively 
small value of the electromagnetic coupling constant en- 
sures that corrections to the one-photon exchange ap- 
proximation should be relatively small. To further im- 
prove on the already impressive accuracy achieved in the 
analysis of electron scattering data, the contribution of 
two-photon exchange (TPE) effects in elastic e-p scatter- 
ing are under theoretical investigation 0, S Hi S IE IE • 
Two-photon exchange effects have garnered particular 
attention as of late due to their potential role in re- 
solving the discrepancy between electromagnetic form 
factors measured through the Rosenbluth Separation 
method | E ITol lll| and a polarization transfer tech- 
nique [IlTill (see Sect.inn. 

Data from elastic and inelastic scattering, both in the 
resonance and deep-inelastic regimes, are here studied in 
light of this concern. There is a newly expanded, sub- 
stantial data set which enables in particular a search for 
non-linearities caused by TPE effects. While dedicated 
measurements have been proposed for elastic data 
and a model-dependent analysis of non-linearity has been 
performed for elastic e-p scattering [Tsj ]. this work rep- 
resents a first detailed, model-independent study of non- 
linear effects in elastic and inelastic e-p scattering data. 



II. ROSENBLUTH SEPARATION TECHNIQUE 

The differential cross section for e-p scattering can be 
expressed in the Born approximation in terms of absorp- 



tion of longitudinal (ctl) and transverse (ctt) virtual pho- 
tons as 



dndE 



7 = F 



aT{W^,Q^) + eaL{W^,Q^) 



(1) 



where is the negative squared mass of the virtual pho- 
ton, is the mass squared of the undetected system, 
and F is the transverse virtual photon flux: 



F = 



aK E' 1 
27r2Q2 ^ 1 _ e 



(2) 



Here, a is the fine structure constant, E and E' are the 
energy of the initial and scattered electron, respectively, 
and K is 



K 



2Mv - Q2 
2M 



(3) 



where M is the mass of the proton and v = E — E' . 
The variable e is the relative longitudinal virtual photon 
flux. Therefore e = corresponds to a purely transverse 
photon polarization. 

The Rosenbluth separation technique is used to sepa- 
rate the longitudinal and transverse components of the 
cross section. Here, Eq. ^ is written in the following 
form: 



1 d2(T 



F d^dE' 



(4) 



In the Born approximation, the left hand side, the re- 
duced cross section, depends linearly on e. To perform 
the Rosenbluth separation, data covering a range in e 
at fixed {W^,Q^) values must be obtained. Any devi- 
ation from linearity must come from higher order terms 
that are not included in the standard radiative correction 
procedures. 
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III. TWO METHODS OF FORM FACTORS 
MEASUREMENTS AND TWO DIFFERENT 
RESULTS 

For the case of elastic scattering, the Rosenbluth sep- 
aration technique is used to extract the form factors Ge 
and Gm, from the e dependence of a reduced elastic cross 
section Ur at fixed Q^, i.e. 



da e(l + t) 

'dn 



^TGlj{Q^) + eGl{Q^) , (5) 



where r = Q^/4:M'^. 

At fixed Q^, the form factors Ge and Gj\/ can be ex- 
tracted from a linear fit in s to the measured reduced 
cross sections. Such a Rosenbluth fit yields tG\,j as the 
intercept and G"^ as the slope. With increasing Q^, the 
cross section is dominated by tG\j, while the relative 
contribution of the G'^ term is diminished. Precise un- 
derstanding of the e-dependence in the radiative correc- 
tions becomes crucial at high values of Q^. Therefore, in 
order to measure the ratio Ge/Gm at high values of Q^, 
a polarization transfer method has also been employed in 
Hall A at Jefferson Lab (JLab). 

In polarized elastic electron-proton scattering, the lon- 
gitudinal and transverse components of the recoil polar- 
ization are sensitive to different combinations of the elec- 
tric and magnetic elastic form factors 0, ■ The ratio 
of the form factors can be directly related to the compo- 
nents of the recoil polarization 



Ge ^ Pt {E + E') tan(ge/2) 
Gm ~ Pi 2M 



(6) 



where Pi and Pt are the longitudinal and transverse com- 
ponents of the final proton polarization, and 9^ is the an- 
gle between the initial and final directions of the lepton. 

Recent measurements from Jefferson Lab using the 
polarization transfer technique to measure the ratio 
Ge/Gm have found that Ge decreases more rapidly than 
Gm at large [13. IT^ . This differs from results ob- 
tained at SLAC in a similar range using the Rosen- 
bluth technique. There exist but two explanations for 
this discrepancy. There is either an unidentified system- 
atic experimental uncertainty in the polarization transfer 
data, or a systematic uncertainty common to all Rosen- 
bluth data. 

It has been estimated that a 5-7% systematic correc- 
tion to the e dependence of the reduced Rosenbluth cross 
section measurements would be needed in order to resolve 
the discrepancy J£,^J. However, a detailed analy- 
sis does not show any inconsistencies in the cross section 
data sets [sj. Moreover , ne w high cross section data 
from Jefferson Lab are consistent with the older 

SLAC cross section data l21l obtained in the same 
range. The results of Ref. [ij, where the struck proton 
rather than the scattered electron was detected, have a 
precision comparable to the polarization transfer mea- 
surements. It must be concluded, then, that the needed 



5-7% e-dependent correction is not due to standard ex- 
perimental considerations in the measured Rosenbluth 
cross sections. 

It has been suggested that the discrepancy may 

be explained by TPE effects not fully accounted for in 
the standard radiative corrections procedure of Mo and 
Tsai [2^ . The polarization transfer technique involves a 
ratio of cross sections, and hence is expected not to be 
very sensitive to such effects 0,0,0- In contrast, these 
contributions can significantly affect the Rosenbluth sep- 
aration technique. TPE contributions can be indepen- 
dent of e (affecting both G^ and in Eq. lO), linear 
in £ (significantly affecting G^), or non- linear in e. 

The experimental evidence for significant TPE contri- 
butions to the form factor measurements is still quite 
limited. While the non-zero transverse beam spin asym- 
metry [2^13 provides direct evidence for TPE in elastic 
e-p scattering, we are lacking similar evidence for such ef- 
fects on the unpolarized cross sections. The discrepancy 
between polarization transfer and Rosenbluth extractions 
of Ge/Gm provides only an indirect indication of a miss- 
ing correction, while direct searches for TPE through the 
comparison of electron-proton and positron-proton scat- 
tering yield some evidence of deviations from the Born 
approximation at low e, but only at the three sigma 
level j25j] . Observation of a deviation from linearity in the 
reduced cross section would provide a clear indication of 
TPE (or other higher order corrections not included in 
standard radiative correction procedures), though only 
the non-linear portion of the correction could be directly 
isolated. New high-precision Rosenbluth data in elas- 
tic [Tl'l and inelastic [2^ e-p scattering allow for a much 
more sensitive search for such non-linearities. 

This work reports results of a search for effects of 
TPE corrections in elastic and inelastic scattering data 
by searching for e-dependent non-linearities in existing 
Rosenbluth separations. We note that this analysis will 
not be sensitive to either systematic shifts in the reduced 
cross section of Eq. Q), or to two-photon effects which 
are linear in e. 



IV. DATA OVERVIEW 

Table HI lists the data sets included in the present 
analysis. We include several measurements of elastic 
e-p to cover a range in , while the SLAC measure- 
ments [23, and the recent JLab measurement 
cover the DIS and resonance regions. 

For elastic e-p scattering, we examine Rosenbluth ex- 
tractions from several different experiments. We study 
the Rosenbluth separation for the experiments and 
values listed in Table II of Ref. 9j , including the updated 
radiative corrections |^ . In addition, data from two re- 
cent Jefferson Lab measurements ^3 are included. 
In all cases, the reduced cross sections are taken from a 
single experiment and single detector. Where necessary, 
cross section values at slightly different values are 
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TABLE I: Summary of experiments included in the analysis, 
including the number of L-T separations and the typical cross 
section uncertainties (excluding normalization uncertainties) . 







#of 


Typ. 


Lab 


Elastic data 


[(GeV/cJ ] 


T T'l. 

L<— is 


QG j G 




Janssens et al. [2^ 


0.2-0.9 


20 


4.7% 


Mark III 


IjlLL CL (jit,. 


Z.(J o.o 


A 


1 7% 
1. 1 /O 




Berger et al. 


0.4-1.8 


8 


2.6% 


Bonn 


Walker et al. [2] " 


1.0-3.0 


4 


1.1% 


SLAG 


Andivahis et al. [21j 


1.8-5.0 


5 


1.3% 


SLAG 


Christy et al. [10] 


0.9-5.2 


7 


1.3% 


JLab 


Qattan et al. LUJ 


2.64-4.1 


3 


0.6% 


JLab 


Inelastic data 


H/2[GeV^] 








Liang et al. [26] 


1.3-3.9 


191 


1.7% 


JLab 


Dasu et al. [27, 28] 


3.2-30 


61 


3.0% 


SLAG 



" Data below 20 deg are excluded. 
* Data from 8 GeV spectrometer. 
Excludes "slope" systematic uncertainties. 
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: W = 2.025 
rQ' = 2.5 

- R = 0.134(63) 


- = 2.175 
-R = 0.349(80) 


- W' = 2.225 
LQ' = 0.75 ^ 

^[^='(1502(9 9) 


- W' = 2.225 
-Q' = 2.5 . 

^1^=^1426(85) 


- = 2.275 
iQ' = 1.75 ,. 

:' R = 0.387(77) 


^W" = 2.325 
LQ' = 0.75 

\ R = 0.175(72) 


Tw' = 2.425 
Lq' = 0.75 
- k^' 

fR = 0.178(70) 


^W' = 2.575 
: = 0.75 

Lr = 0.218(70) 


^W' = 2.625 
: Q' = 0.75 

--R = 0.226(64) 


; = 2.675 
: Q' = 0.75 

i R = 0.252(63) 


; W' = 2.725 
;Q' = 2.5 

r 

Ar = 0.353(74) 


W' = 2.775 
^^Q^j^Z5^, 

- R = 0.402(76) 

,,,1,2,,' 


IW' = 2.875 
- = 1.75 

"'r^O. 348(55) 


: W' = 2.925 
: = 0.75 

AR = 0,222(57) 


rW = 2.975 
LQ' = 0.75 

Tr = 0.222(56) 


-W' = 3.05 
: = 0.75 

Ar = 0.316(86) 



0.5 0.5 0.5 0.5 



shifted to a fixed values. Only small corrections were 
needed, typically below 2%, although a handful of points 
were corrected by 5-10%. There are a total of 51 Rosen- 
bluth separations that we will examine for non-linearities. 
Typical point-to-point uncertainties are roughly 1-2% for 
most of the data sets, although several of the older experi- 
ment had larger uncertainties and the EOl-OOl results [Tlj 
have point-to-point uncertainties below one percent. 

For the resonance region, we used newly obtained data 
from JLab Hall C experiment E94-110 26, 32], which 
was utilized to separate the longitudinal and transverse 
unpolarized proton structure functions in the nucleon res- 
onance region via the Rosenbluth separation technique. 
The experiment ran with seven different energies rang- 
ing from 1.2 GeV to 5.5 GeV over a scattering angle 
range 12.9<6'e<79.9. The total point-to-point uncer- 
tainty on the cross section measurements was approxi- 
mately 2% [s^- The data taken from this experiment 
were used to perform 191 Rosenbluth separations cov- 
ering the kinematic region 0.5<Q^<5.0 (GeV/c)^ and 
1.1<W^<4.0 GeV^. Examples of these Rosenbluth sepa- 
rations are shown in Fig. ^ These data were used to ex- 
tract the ratio, i?, of longitudinal to transverse cross sec- 
tion components. Rosenbluth separations are performed 
in five bins and 43 bins. The cross section val- 
ues are interpolated to the central W"^ and values of 
each bin u sing a global fit to the world's resonance re- 
gion data |26l l33l|. with constraints built in to provide 
a smooth transition to the DIS region and the ^ 
limit. Typical corrections were 5%, and data that re- 
quired corrections larger than 60% were excluded. While 
a few points had large enough corrections that the uncer- 
tainty in the interpolation may not be negligible for the 
given data set, their effect on the extracted values of P2 
should be small and random, thus providing a negligible 
contribution to the uncertainty in the extracted limits on 
P2. 



FIG. 1: (Golor online) Example Rosenbluth separations per- 
formed in experiment E94-110. Each figure includes the kine- 
matics, W"^ in GeV^ and in (GeV/c)^, and the extracted 
value of R = gl/ o"t . 



rW"=6.55 
rQ'=0.63 ,j 

HR=0. 324(94) 


rw'=9.88 

rR=0, 382(87) 


hW'=12.66 
|:r=0. 344(96) 


hW'=15.88 
^Q^=5 

rR=0. 173(69) 


lW"=9.19 
: Q'=4 

^R=0, 174(70) 


UW'=4.15 
^Q"=4 
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W'=9.06 

: Q'=10 

^K'-* 

LR=0. 093(65) 


LW'=10.48 
lQ"=16 

[r=0. 278(1 12) 
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FIG. 2: (Golor online) Example Rosenbluth separations per- 
formed in experiment E140 at SLAG |2||. Each figure 
includes the kinematics, IF^ in GeV^ and in (GeV/c)^, 
and the extracted value of i? = (jl/ctt. 



For the deep-inelastic scattering (DIS) region, the 
data from experiment E140 at SLAG [23, 123 have 
been used. A total of 61 Rosenbluth separations 
have been performed covering the kinematic region 
0.63<Q2<20 (GeV/c)2 and 2.5<W2<30 GeV^. The to- 
tal point-to-point uncertainties on the cross section mea- 
surements depends on the actual kinematics, but are typ- 
ically 2-3%. Example Rosenbluth separations from E140 
are shown in Fig. [2 

The data from the JLab E94-110 and SLAG E140 ex- 
periments are by far the largest data sets available for 
inelastic Rosenbluth separations, while the new elastic 
measurements from JLab EOl-OOl provide significantly 
more precise data for elastic scattering. 
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V. ANALYSIS AND RESULTS 

From the discussions of Sections and IIIII it is clear 
that the Hnearity of the reduced cross section in Eq. Q 
is a crucial component of the Rosenbluth technique, and 
that two-photon exchange corrections could introduce a 
non-linear e-dependence in Eq. (0J. Therefore, such cor- 
rections could manifest themselves as non-linearities in 
Figs. [2 and El 

To search for such non-linearities, the following anal- 
ysis has been performed. For each data set with three 
or more e values at fixed and VF^, the reduced cross 
sections are fit to a quadratic in e, of the form 

(7, = Po ■ [1 + Px{e - 0.5) + P2(e - 0.5)2]. 

In the absence of TPE, we expect Pq = ctt + 0.5(Tl, 
Pi = c7l, and P2 = 0. TPE corrections can modify Pq 
and Pi, and may introduce a non-zero value of P2, the 
fractional curvature relative to the Pq, the cross section 
at e = 0.5. The only estimates we have for the size of 
the non-linearity come from calculations for elastic e-p 
scattering. If one takes the calculations 0) 13 TPE 
effects for elastic scattering and scales the size of the 
corrections so that they are large enough to explain the 
discrepancy betwe en p olarization and Rosenbluth extrac- 
tions, as done in one obtains P2 values of «6-9%, 
although the precise value depends significantly on 
and the e range of the data. 

While P2 represents the fractional curvature, the size 
of cross section deviations from linearity will be much 
smaller. For P2 = 10%, the maximum deviation of the 
cross section from P2 = would be 2.5%, at e = 0,1. 
The effects are even smaller if the e range of the data, 
Ae, is less than one. In this case, the deviations from 
P2 = will go approximately as (Ae)^. Finally, when 
one performs the Rosenbluth separation, the extracted 
values of and ot will be shifted from their true values 
in order to minimize the deviation from the straight line 
fit, reducing the deviations by roughly a factor of two 
from the size of the P2 contribution in Eq. [7| Thus, the 
maximum observed deviations from linearity will be, 

A_. = ^^-^i^^raa. ^ ^^^^2/8^ 

a 

typically more than a factor of ten smaller than the value 
of P2. For the expected P2 values of <^10% and a rather 
large Ae range of 0.8, one expects maximum deviations 
from linearity to be at the level of <0.8%. So even with 
high precision measurements and a good e range, one 
needs a large data set to provide meaningful limits on 
the non-linearities. 

We perform the fit from Eq. [7|for each of the elastic, 
resonance region, and DIS Rosenbluth data sets. Fig- 
ures 13 and 0] show P2, binned in for the elastic data, 
and binned in W for the resonance and DIS data. The 
results are consistent with no non-linearities, and there 
is no apparent or dependence. Table UTI shows the 
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FIG. 3: (Color online) Extracted values of the curvature pa- 
rameter, P2, as extracted from the elastic data as a func- 
tion of . The red dotted line indicates the average value, 
(P2> =0.019 ±0.027. 
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FIG. 4: (Color online) Extracted values of the curvature pa- 
rameter, P2 , as extracted from the inelastic data as a function 
of W . Data in each W bin is averaged over all values in 
the resonance region and DIS measurements. The red dotted 
line indicates the average value, (P2) = —0.048 ± 0.036. 



extracted value for P2, the 95% confidence level upper 
limit on I P2 | , and the approximate maximum deviation 
from linearity for the elastic, resonance region, and DIS 
{W^ > 4 GeV^) data. From these results, we determine 
the 95% confidence level upper limits on | P2 | 0/6. 4% 
for the elastic data and 10.7% for the inelastic data. This 
yields limits on the deviations of the data from the Rosen- 
bluth fit of roughly 0.4% (0.7%) for the elastic (inelastic), 
assuming a Ae range of 0.7. 

Note that it is also possible for a purely linear correc- 
tion to introduce a small non-linearity, since 

— = (1 + Pe) • (1 + C2-^e) = 1 + (P + C2-y)£ + R ■ C2^e^ 

(9) 

where P — (Jl/ctt and (1 -I- C2-yS) is the TPE correc- 
tion factor. For the elastic data at high and all of 
the inelastic data presented here, P <^ 0.2 — 0.3, while 
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TABLE II: Extracted values and 95% confidence level upper 
limit on B2. A,„ai is the upper limit on deviations of the 
cross section from linearity (Eq.jSJ. 







1 P2 \max 

95% C.L. 


^max 

95% C.L. 


Elastic 

Resonance 

DIS 


0.019(27) 
-0.060(42) 
-0.012(71) 


0.064 
0.086 
0.146 


0.8%-(Ae)^ 
l.l%-(Ae)2 
1.8%-(Ae)2 



estimates TPE predict a change in slope, 6*2^, of approx- 
imately 5%. Hence, the nonlinear term arising from a 
Unear TPE correction, R ■ will be very small, yield- 
ing P2 ^ 1%. At low values, the value of R for the 
elastic cross section becomes quite large, yielding val- 
ues of P2 on the scale of C2-y for R ^ 1. However, R 
is only this large for ^ o.4 GeV^, where the TPE 
corrections decrease as ^ |E [H El . The ef- 
fect is ^1% if one assumes that C27 increases slowly as 
one goes up from = 0, as it does in calculations 
and phenomenological extractions of the TPE correc- 
tions 19, ^J. If one takes a more rapid increase with 
Q^, C2-y = 0.06- [1 -exp(QV0.5GeV2)], we obtain values 
of P2 coming from the linear correction of 1.5-2.5% for 
< I GeV^. Thus, the size of this effect should be well 
below the sensitivity of the existing data in all cases. 

To better visualize the limits on nonlinearities, we have 
also performed a global comparison of the residuals be- 
tween the reduced cross sections and a linear fit to the 
reduced cross sections. For the previous fit, data sets 
with a very small Ae range have very little sensitivity to 
the curvature. Although these data sets have large uncer- 
tainties, they still provide meaningful P2 values. When 
plotting the residuals, we want to exclude such data sets 
because the data points have uncertainties comparable 
to the other data sets, but the residuals little sensitivity 
to non-linearities. Thus, we require include only those 
data sets where Ae > 0.4 when studying the residuals. 
This cut reduces the number of data sets to 35 for elastic 
kinematics, 77 in the resonance region, and 38 in the DIS 
region. 

For the data sets with sufficient e coverage, we take 
to be the residual from the Roscnbluth (one-photon 
exchange) fit. 



— 



O'Data — <y fit 



(10) 



and obtain a value of Ri~^ for every cross section mea- 
surement in the Rosenbluth data sets. We can then de- 
termine the weighted average value from the extracted 
Rij values in e bins for the elastic, resonance region, and 
inelastic data sets. In the absence of TPE contributions, 
one expects Ri^ = in every e bin and, hence, any s 
dependence to Ri-y is an indication of two-photon ex- 
change. Figure [5I shows the combined Rij as a function 
of e for the elastic, resonance region, and DIS data, and 
all three data sets are consistent with Rij=0. One can 
see that while the elastic and resonance region data have 
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FIG. 5: (Color online) The weighted average of Ri-y = 
[oData — Fit) I a Fit for the elastic measurements, the reso- 
nance region data (JLab experiment E94-110) and the DIS 
data (SLAC experiment E140). The solid lines are the fits to 
the form of Eq. [TTl 



high precision, the resonance region data has less data as 
£ and 1 , while the DIS data has limited data an low 
£, as well as lower statistical precision in general. 
We fit the combined residuals to the form 



i?!-, = A + B(£ - £0)^ 



(11) 



where B, and £0 are the fit parameters. Because we 
are fitting to residuals that have already had the expected 
linear cross section dependence removed, we expect that 
R\^ will average to zero, yielding A w in the absence of 
any strong £ dependence. Indeed, we find A <, 0.05% for 
the elastic, resonance, and DIS data. The quadratic fit 
to Ri~^ yields a curvature parameter, B, consistent with 
zero for all data sets. We obtain B = (0.9 ± 2.0)% for 
the elastic data, (—2.3 ± 3.0)% for the resonance region 
data, and (0.9 ± 3.8)% for the DIS measurements. 

While the limits in Tab. ^1 provide the best quanti- 
tative limits on deviations from linearity, the residuals 
shown in Fig. Ogive a better idea of the sensitivity of the 
different data sets in different regions of s. The parame- 
terization of Eq.EI assumes a simple quadratic non-linear 
term, while some models for the contribution to elastic 
scattering indicate larger non-linearities for £ — > 1. From 
Fig. [51 we see that this region is not as well constrained 
for the resonance region data, while very low e values are 
not well constrained except in the elastic data. 
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VI. CONCLUSION 

We have searched for possible two-photon exchange 
contributions that show up as non-hncaritics in Roscn- 
bluth separations. We have used existing data in the elas- 
tic and deep-inelastic scattering region and recent data 
in the nucleon resonance region. Wc do not find any ev- 
idence for TPE effects. The 95% confidence level upper 
limit on the curvature parameter, P2, was found to be 
6.4% (10.6%) for the elastic (inelastic) data. This limits 
maximum deviations from a linear fit to ^0A% (0.7%) 
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